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We discuss the possibility that r-process nucleosynthesis may occur in the winds from 
gamma ray burst accretion disks. This can happen if the temperature of the disk is 
sufficiently high that electron antineutrinos are trapped as well as neutrinos. This implies 
accretion disks with greater than a solar mass per second accretion rate, although lower 
accretion rates with higher black hole spin parameters may provide viable environments 
as well. Additionally, the outflow from the disk must either have relatively low entropy 
s ~ 10 or the initial acceleration of the wind must be slow enough that it is neutrino and 
antineutrino capture as opposed to electron and positron capture that sets the electron 
H ' fraction. 



1. INTRODUCTION 



X ■ The r-process of nucleosynthesis is responsible for over half the elements with A > 100 
but to date no self-consistent model for an astrophysical site for the production of these 
elements has been identified. The two candidates most discussed in the literature are the 
neutrino driven wind of Type II supernovae jl!2j and ejecta from neutron star mergers 
The former site is attractive because it occurs on a timescale such that it could 



reasonably account for observations in metal poor halo stars e.g. [S| and could produce 
roughly the right amount of material [2|. Furthermore, the neutrinos from the protoneu- 
tron star set the electron fraction to a relatively low value and this, combined with the 
high entropy, makes it a promising site. Although attempts to produce the r-process 
elements this way initially succeeded, subsequent refining of the models indicated that 
the conditions were a near miss, e.g. |7|8j . In addition, self-consistent inclusion of the 
neutrinos in a reaction network demonstrated that electron neutrino capture on neutrons 
during the course of alpha particle formation drove the electron fraction up to values 
unacceptable to the r-process 0. 

The other primary candidate for the production of r-process elements, the ejecta from 
neutron star mergers, is a viable candidate for producing some of these elements. However 
as a dominant r-process site, it would be difficult to reconcile with the observations from 
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metal poor stars, e.g. [TU]. If the r process has more than one component, as suggested by 
meteoritic evidence [llj and observations of the low mass end of the r-process distribution, 
then neutron star mergers could well contribute. 

Here we discuss another possible site for the r process, which is in the winds from 
accretion disks in gamma ray bursts. Gamma ray bursts represent an emerging new 
area in nucleosynthesis research. The leading candidate sites for these objects are rare 
supernovae, such as a collapsars ^2], and for the short bursts, neutron star-neutron star 
mergers |13|14j . In either case, the system likely forms an accretion disk surrounding a 
black hole. Some nucleosynthesis will occur in the jet [T^], some will occur as explosive 
burning (TE] and still more will occur in the outflow, perhaps a wind |12|17j . from the 
accretion disk. The rate of gamma ray bursts is estimated to be quite low, ~ 1CT 5 
per year in the Galaxy, e. g. as discussed in However since their origin is not 

yet fully understood, they may well represent one end of a continuum of objects. It is 
therefore difficult to estimate the rate of ejection of their nucleosynthesis products into 
the interstellar medium. 




Figure 1. Electron neutrino and antineutrino surfaces in a M = 1 M / s, a = (a, on the 
left) and a M = 1OM / s, a = (b, on the right) accretion disk. The solid line shows the 
density scale height of the disk. 



The possibility of obtaining an r process in an accretion disk wind is closely tied to 
understanding the neutrinos which are emitted copiously from the disk. Neutrinos are 
involved in all charge- changing interactions 

e~ + p <-> v e + n (1) 

e + + n <-> v e + p (2) 

both in the disk and in the outflow from the disk. In order to make the outflow neutron 
rich, either electron capture or electron antineutrino capture must dominate. In the case 
of disks with high accretion rates where both the neutrinos and the antineutrinos are 



3 



trapped (as in Fig. [IJ, it is antineutrino capture that can create conditions conducive to 
forming the r-process elements. 

2. NEUTRINO TRAPPING IN THE DISK 

Whether the neutrinos in the accretion disk are trapped or not depends on the param- 
eters that characterize the disk, such as the accretion rate M, the viscosity a and the 
black hole spin parameter a. In Fig. Q we show the trapped regions for neutrinos and 
antineutrinos which we have calculated from the disk models of DiMatteo, Perna and 
Narayan (2002) |19j . In the outer regions of the disk, neutrinos and antineutrinos are 
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Figure 2. Electron neutrino (solid line) and antineutrino (dashed line) spectra at a point 
well above the decoupling surfaces, z = 220 km and r = 250 km in a M = 1 M / s, a = 
(a, on the left) and a M = 10 M Q / s, a = (b, on the right) accretion disk. Note that the 
y-axis scale is different in the left and right panel. 

produced by the inverse beta decay processes (forward reactions in Eqs. ^andEJ) and the 
neutrinos escape freely. However, toward the center of the disk, the electron neutrinos 
first become trapped, and then the electron antineutrinos become trapped as well, en- 
abling the backward reactions in Eqs. ^and|2]to come into equilibrium with the forward 
processes. In the higher accretion rate model, M = 10 M / s, a = these regions are 
larger than for the more moderate accretion rate model M — 1 M Q / s, a = 0. In particular 
the region enclosed by the antineutrino surface has expanded, although the antineutrinos 
still decouple considerably further in than the neutrinos. More detail on calculations of 
the disk neutrinos can be found in |20j . 

We use the neutrino and antineutrino spectra emitted from every radius on the disk 
surface to determine the neutrino and antineutrino fluxes at every point above the disk. 
The spectra from the electron neutrinos and antineutrinos at about 200 km above the two 
disks shown in Figs, and^o can be seen in Figs. |2k and|2fc>. For the lower accretion rate 
disk, the electron neutrino flux far dominates the electron antineutrino flux. Although the 
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electron antineutrinos which come from the surface shown in have higher temperature 
than the neutrinos, the area of the antineutrino surface is very small, so the number 
of antineutrinos is not large and a significant number of those seen in |2k are coming 
directly from inverse beta decay in the free streaming region. In the higher accretion rate 
disk, the situation is quite different. Because the antineutrino surface has grown, and 
the antineutrinos have higher temperature than the neutrinos, the antineutrino flux is 
larger and more energetic than the neutrino flux. In these models the electron neutrino 
temperature at the neutrino surface varies with position but is around T MeV ~ 2.5 to 
TmsV ~ 4.5 while the temperature at the antineutrino surface is around TueV ~ 3.6 to 
TMeV ~ 5.1 |2*U] . In even lower accretion rate disks (not shown here), e.g. Popham, 
Woosley, Fryer (1999) [2T|, M = 0.1M Q / s, a = 0.95, the neutrinos are barely trapped 
and the antineutrinos are not trapped at all. 

3. OUTFLOW FROM THE DISK 

We study the effect of the charge-changing interactions on the outflow from the disk by 
using a parameterization for the wind, as in |22] • The velocity is taken to be 

\u\ = Voo (l - (3) 

where R = (z 2 + r^) a5 for the first, vertical part of the trajectory and the starting position 
of the material is Rq. Here z is the vertical coordinate above the disk, and r c is the 
cylindrically radial component along the disk. The parameter f3 controls the acceleration 
of the wind; for lower f3 the wind accelerates faster. In terms of the neutrinos, larger (3 
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Figure 3. Electron Fractions as a function of distance for two different disk models M = 
1 M / s, a = (a, on the left) and M = 10 M / s, a = (b, on the right) are shown. The 
conditions in each case are s = 40, r$ = 250 km, (3 = 2.5, Voo = 3 x 10 4 kms _1 (dashed 
line) and s = 40, r = 250km, j3 = 0.8, v oo = 3 x 10 4 kms~ 1 (solid line). The dotted and 
dot-dashed lines show the same calculation without the neutrinos included. 
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means more time for the neutrinos to influence the composition of the outflow. In Fig 
El we show outflow from two different types of disk models, DPN M = 10 M Q / s and 
M = lMg/s. 

It can be seen from these figures that the neutrinos influence the electron fraction 
in opposite directions in the two models. In the more moderate accretion rate model 
M — 1 M Q / s, a = 0, the large flux of neutrinos shown in Fig praises the electron fraction 
to quite high values. However in the higher accretion rate model, M = 10 M / s, a = 0, 
the large flux of antineutrinos decreases the electron fraction to quite low values, even in 
the case of a high entropy outflow. Also interesting is the low entropy case in Fig. Eb, 
where the electron fraction remains quite low regardless of the influence of the neutrinos. 
This is due to the electron degeneracy of the material at a high density and low entropy. 

In Fig. HJ we survey a number of different outflows, by plotting the electron fraction 
measured at a little under an MeV against entropy for several different values of (3. In all 
cases the final outflow velocity is = 10 4 kms _1 . 




Figure 4. Electron fractions in the outflow measured at a temperature of 1 MeV for two 
different disk models M = 1M / s, a = (a, on the left) and M = 10M Q / s, a = (b, on 
the right) are shown. Various values of the acceleration parameter j3 are plotted against 
entropy per baryon in the outflow. 



The plus signs show the effect of a slow outflow parameter (3 — 2.5. In this case the 
neutrino and antineutrino capture rates completely overwhelm the electron and positron 
capture rates, and the system finds a weak equilibrium such that 

Y e = (4) 

where A^ eP is antineutrino capture on protons and X Uen is neutrino capture on neutrons. 

Again we see that a slow outflow from a high accretion rate disk will produce a low 
electron fraction due to the neutrino interactions, although for any outflow a low entropy 
will also create a low electron fraction. For the lower accretion rate disk and slow outflows, 
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the neutrinos always make the material more proton rich and in fact it is only the case of 
low entropy and minimal neutrino interactions that produces a low electron fraction. 



4. NUCLEOSYNTHESIS 

Although there are two regions of parameter space which will produce very neutron rich 
winds, this is not a guarantee of a successful r process, so in this section we investigate 
this possibility. We take the high accretion rate model M = 10 M Q / s, a = and an 
outflow with a fairly low entropy, s=10 and (3 = 0.8 and a final velocity of 3 x 10 4 kms _1 . 
This is the same trajectory as shown as the solid line in Fig. 03. Our calculation indicates 
that the low entropy trajectories do produce an r process as shown in the solid line in 
Fig. EJ The dashed line shows a calculation with a higher entropy, s = 40, but a slower 




Figure 5. Shows two r-process nucleosynthesis calculations, both in outflows from a 
M — 10 M / s, a = accretion disk. We show a low entropy s = 10 outflow with (3 = 0.8 
and Voo = 3 x 10 4 kms _1 as the solid line and we show a high entropy outflow, s = 40 but 
with a slow acceleration parameter, (3 = 2.5 and = 3 x 10 4 kms _1 as the dashed line. 
The crosses show the measured solar abundances. 



7 



outflow acceleration (3 = 2.5 so that the neutrinos have greater influence. Note that in 
neither of these calculations is there an "alpha effect" |8|23| and this is due to the high 
outflow velocity at the time alpha particles are forming. 

While the higher accretion rate disks may be a viable site for the r process, moderate 
accretion rates disks such as M = 1M Q / s, a = will produce much higher electron 
fractions. In Fig. 0^, it can be seen that the material which has been irradiated by a 
large neutrino fluence because of its slow acceleration ((3 = 2.5) is very proton rich. This 
type of outflow will produce primarily nickel, but it will also produce some nuclei on the 
proton rich side of the valley of beta stability as well, such as 58 Cu, 59 Zn, 50 Fe and 52 Fe. 
For still lower accretion rate models, such as M = 0.1 M Q / s, a = 0.95 the outflow will 
have an electron fraction which is closer to Y e = 0.5, although the neutrinos may drive the 
electron fraction up to as much as Y e = 0.6. Winds of this type with entropies of order 
s ~ 30 may have an unusual nucleosynthesis pattern with large overproduction factors of 
elements such as 42 Ca and 45 Sc, 46 Ti, 49 Ti, 63 Cu, 64 Zn as discussed in Lower spin 
parameter models are discussed in |2"5] . 

5. CONCLUSIONS 

The winds from accretion disks surrounding black holes in the context of gamma ray 
bursts are a new arena in which to investigate nucleosynthesis. The most important 
parameters which determine the elements formed in the outflow from the disk are the 
accretion rate and black hole spin parameter. This is because the density and temperature 
of the disk determine where the neutrinos and antineutrinos become trapped. Neutrinos 
become trapped before antineutrinos, and in some disks only neutrinos are trapped, not 
antineutrinos. In these disks the wind is proton rich, from the reaction u e +n — *> e~+p, and 
produces considerable nickel-56, but also elements on the proton rich side of the valley 
of beta stability. However as the accretion rate increases, a sizable region of trapped 
antineutrinos can develop, and these antineutrinos cause the wind to become neutron 
rich. This happens because the higher temperature antineutrinos cause v e + p — > e + + n 
to be faster than the corresponding reaction for neutrinos. If the outflow is is slow enough, 
for example in our parameterization at around (3 = 2.5, or has low entropy s ~ 10, then 
the electron fraction is sufficiently low that r-process elements are produced in the outflow. 
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